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Abstract 

The Riesz transform is a natural multi-dimensional extension of the 
Hilbert transform, and it has been the object of study for many years 
due to its nice mathematical properties. More recently, the Riesz trans- 
form and its variants have been used to construct complex wavelets and 
steerable wavelet frames in higher dimensions. The flip side of this ap- 
proach, however, is that the Riesz transform of a wavelet often has 
slow decay. One can nevertheless overcome this problem by requiring 
the original wavelet to have sufficient smoothness, decay, and vanish- 
ing moments. In this paper, we derive necessary conditions in terms 
of these three properties that guarantee the decay of the Riesz trans- 
form and its variants, and as an application, we show how the decay 
of the popular Simoncelli wavelets can be improved by appropriately 
modifying their Fourier transforms. By applying the Riesz transform 
to these new wavelets, we obtain steerable frames with rapid decay. 

1 Introduction 

The Riesz transform of a real- valued function f{x) on R'^ is the vector- valued 
function Uf{x) = (7^l/(a;), . . . ,Tldf{x)) given by 

n,f{x) = Ca lim / /(y)^i^^dy, (1) 

^^0 7||a:-y||>e \\x - yY+^ 
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where Cd is a absolute constant depending on the dimension d [5j. In this 
formula, we can see that the Riesz transform is defined by a principal value 
integral, as is the case for the Hilbert transform. Roughly speaking, the i-th 
component of the Riesz transform is obtained through the convolution of 
f{x) with the kernel Xj/||a;||'^"'"^, where ||a;|| is the Euclidean norm of x, and 
Xi are its coordinates. The truncation of the integral in ([1]) is used to contain 
the singularity of the kernel at the origin, by systematically "chopping" it 
off around the origin. This added detail makes the analysis of the transform 
quite involved. 

Though the transform appears arcane at first sight, it tends to emerge 
naturally in various physical and mathematical settings due to its funda- 
mental nature. It is essentially the unique scalar-to-vector valued transform 
that is unitary (conserves energy), and is simultaneously invariant to the 
fundamental transformations of translation, scaling, and rotation. For ex- 
ample, the role of the Riesz transform in the theory of differential equations 
and multidimensional Fourier analysis has long been recognized in mathe- 
matics, where the above-mentioned invariances play a central role, see, e.g., 
[51 [20| I21j . The invariance of the transform to translations, namely that 
TZ\f{- — u)][x) = [JZf]{x — u), is immediately seen from ([T|). On the other 
hand, its invariance to scaling is confirmed by the change of variable x sx 
in ([1]), where s > 0. The other properties are, however, not obvious at first 
sight. Instead, these are better seen from its Fourier transform, which is 
given by (see, e.g., fSOj) 



From this and Parseval's theorem, we see the sum of the squared L^-norms 
of TZif{x) equals the squared L^-norm of f{x). That is, the energy of f{x) 
is divided between the different Riesz components. The Fourier expression 
also tells us that f{x) can be recovered from its Riesz transform using a 
simple reconstruction, namely as 



Using the invariance of the Fourier transform with rotations, it can also 
be verified from ([2]) that the Riesz transform of the rotation of f{x) is the 
(vector-field) rotation oiTZf{x). 

The Riesz transform is a natural multidimensional extension of the Hilbert 
transform, obtained by setting d = 1 in ([1]). The fact that the Hilbert trans- 
form plays a central role in determining the (complex) analytic counterpart 
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f{x) = -Y,T^i[T^if{^)]- 



i=l 
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of real-valued signals, and, particularly, in determining its instantaneous 
phase and frequency, has long been known in the signal processing and 
wavelet community, e.g., see [U HI O [T7|. The introduction of the Riesz 
transform to the community, however, has been more recent. In their pi- 
oneering work m [5] , Felsberg and Sommer used the transform to define a 
two dimensional generalization of the analytic signal, which they called the 
monogenic signal. At around this time, Larkin et al. introduced a complexi- 
fied version of the Riesz transform in two dimensions [141 [T5] . They showed 
that this complex transform could be used for demodulating interferrograms, 
and for analyzing fringe patters in optics. Felsberg's work on the monogenic 
signal triggered a series of research works relating to the applications of the 
Riesz transform in the design of wavelets. For example, this inspired Metikas 
and Olhede to realize a monogenic version of the continuous wavelet trans- 
form in [16]. Further work, concerning filters based on Riesz transforms, can 
be found in [H [El |13] . 

More recently, in a series of papers by Unser et al. [22l [231 121]) it was 
shown how certain functional properties of the transform, particularly its 
invariances to translations, scaling and rotations, and its connection with 
derivatives and the Laplacian, could be used for realizing complex (mono- 
genic) wavelet bases and steerable wavelet frames in higher dimensions. This 
provided a natural functional counterpart to Simoncelli's steerable filter- 
banks |19j . which is a popular framework for decomposing images at mul- 
tiple scales and orientations. The key connection provided by Unser et al. 
was that the steerable wavelets of Simoncelli could be realized as the Riesz 
transform of appropriate bandpass functions. This powerful functional inter- 
pretation allowed them to enrich the wavelet design along several directions. 
For example, the authors proposed the so-called higher-order Riesz trans- 
forms obtained through repeated applications of the various components of 
the Riesz transforms. They showed these higher-order transforms could be 
applied to a given wavelet to obtain a larger palette of steerable wavelets. 
When applied on the data, these frames act as (multi-scale) partial deriva- 
tives. The important point was that all of this could be done using fast and 
stable filterbank algorithms. 

An aspect of this theory that has not previously been studied in depth 
is how the Riesz transform affects the decay of a wavelet. The importance 
of this issue lies in the fact that a slowly decaying wavelet can produce 
poor approximations that suffer from problems such as ringing artifacts and 
boundary effects. Here, we propose that good decay of the Riesz transform 
TZf{x) is determined by the smoothness, decay, and vanishing moments of 
f{x). In fact, this observation explains the poor spatial decay of some of 
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the earlier wavelet constructions involving the Riesz transform. For example, 
Felsberg applied the transform on the Poisson kernel which has no vanishing 
moments [7[[8], and Simoncelli used it on wavelets with poor spatial decay 
(due to the non-smooth cut-off in the spectrum) [19j . 

In this paper, we perform a quantitative study of the decay property of 
the Riesz transform and its higher-order variants. We expand upon the work 
in [3], where the intimate connection between the Hilbert transform (the one- 
dimensional Riesz transform) and wavelets was investigated. It was observed 
there that for a wavelet with large number of vanishing moments and a 
reasonably good decay, its Hilbert transform automatically has an ultra-fast 
decay. We extend the results in [3j first to the Riesz transform, and then to 
its higher-order counterparts. While the main principles are the same, we are 
required to use a different approach in the higher dimensional setting. The 
reason for this is that when d = 1, the Riesz multiplier has the particularly 
simple form l/vrx, which does not involve the term Also, the concept 
of higher-order transforms is trivial for d = 1, since by applying the Hilbert 
transform twice, we get back the negative of the original wavelet. The higher- 
order transforms are thus of significance only for d > 2. The analysis of these 
transforms presents additional technical challenges due to the asymmetry 
between the decay of the original wavelet and its Riesz transform. As a 
consequence, the results for the higher-order transforms cannot simply be 
obtained by a direct application of the results for a single Riesz transform. 

2 Notation 

We will follow the notations in \22\ [23} . The components of the higher- 
order Riesz transforms are obtained by iterated applications of the compo- 
nent transforms. Thus, given a multi-index of non-negative integers a = 
(«!, . . . , ad), we define the higher-order Riesz transform by 



where T^f denotes the A:-fold application of TZi to f{x), and |q| is the order 
of the transform. Here we are using the the notation \a\ = ai -\- ■ ■ ■ + 
and a! = ail ■ ■ ■ a^!. Also, we will write the partial derivative D"^ f{x) as 




f)°'^ . . . fl"'^ 
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for a given multi-index a = {ai , ■ ■ ■ , aa) ■ The inner-product between two 
vectors x and y will be denoted hy x-y, and a;" will represent the monomial 

It is clear that the various components of TZf{x) have similar analytic 
properties. We therefore fix i once and for all and define the Riesz kernel 
R{y) = 2/i/l|y ll'^"''^- Out goal in this paper is to study the Riesz transforms 
of well-behaved functions, with a certain number of vanishing moments. By 
well-behaved, we mean that they should be sufficiently smooth and must 
have enough decay. To this end, we define the class Ljv(R'^) to be the 
collection of differentiable functions f{x) such that 

1. \f{x)\<C{l + \\x\\)-''-^+^, 

2. \D'^f{x)\<C{l + \\x\\)-'^-^-^+', \a\ = l, and 

3. Jxl^f{x)dx = 0, \P\ < N. 

Here C > and < e < 1 are arbitrary constants, and the class index AT > 1 
is an integer. The first two properties address decay and smoothness, while 
the third condition guarantees vanishing moments. 



3 Riesz transform of wavelets in L]\j(Il^) 

This section comprises the heart of the paper, as we establish results showing 
how the component Riesz transforms preserve decay and vanishing moments 
for functions of the above class. In particular, we will show that the trans- 
form TZif of a function / G LAr(R"') has a rate of decay comparable to / and 
has the same number of vanishing moments. Significantly, these theorems 
are formulated in a way that allows then to be applied iteratively so that 
we may derive results for the higher-order transforms in the next section. 
However, as the transform is given in terms of a singular integral, before 
estimating its decay, we need to show that it is well-defined. 

Proposition 3.1. For all f G -Ljv(R'^), 7^^/ is well-defined, continuous, 
and bounded by 

cfsupiiv/(t)ii+ii/ii,y 

for some C > 0. 
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Proof. We can see that the domain of the singular integral defining TZif is 
unbounded, and wc have an additional difficulty in a neighborhood of the 
point X. We shall apply different techniques to handle these problems, so 
we split the integral into two terms. For a fixed x, we have 



j\\x-y\\>\ II"'' y\\ 

Using the fact that the Riesz kernel is odd and the inequality 



\f{y)- f{x)\\Xi-yi\ 1 1 

||a; — y|| ||a; — y|| ||a; — y||" ^ ^^-^d \\x — y\\°- ^ 

holds, we can write 



Jl>||a;-y||>e 11^ VW Jl>||a;-y||>e 11"^ UW 

(4) 

and apply the dominated convergence theorem, which shows that the first 
term is convergent. This also shows that the first term is bounded by 
Csuptgjtd||V/(t)||. For the second term, we have 



\x — yll'^"'"''' 



which implies it is bounded by ||/||;^. 

To finish the proof, we must show that TZ-if is continuous, so we fix 
X G R"^ and 7 > 0. Continuity can then be proved by showing that there 
are finite constants C > and (5 > such that \1Zif{x) — lZif{z)\ < C7 for 
II a; — z II <5. To that end, consider the transform of / written as 



■if{x) = lini / f{y)- — - — rArrdy 
+ I /(^) || TLi dy 

j\\x-y\\>t' y\\ 
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for e' > e > 0. Using the same techniques as in equations ^ and we can 
see that e' should be chosen so that the absolute value of the first integral 
is smaller than 7. Therefore 



\n,f{x)-nif{z)\<2j+ I \f{y)\ 

\\x—y\\>t' 



L 



Zi - Vi 



\x-y 



\d+l 



\z-y 



\d+l 



and since the Riesz kernel is uniformly continuou^ away from its singularity, 
there is some 6 > (depending on e') for which 



\f{y)\dy 



x—y\\>e' 

<7(2+||/||i) 



when \\x — z\\ < 5. □ 

Now that we have established the fact that the transform is well-defined, 
our focus is shifted to the main result concerning decay. Our approach is 
similar to the one used by Adams in [2j in order to bound weighted norms 
of Riesz transforms, and as in that paper, a key property is the number of 
vanishing moments. As there are different issues with the singular integral 
depending on the region of R*^, we divide it into three pieces and bound 
each separately. One of these regions will be handled by using the vanishing 
moment condition on / and approximating the shifted version of the Riesz 
kernel R{x — •) by its Taylor expansion at 0. Therefore, for a fixed x, we 
denote the degree — 1 Taylor polynomial approximant to R{x — •) by 
Pn{x,-) and define ^]\f{x,y) := R{x — y) — PN{x,y). With this notation, 
we have 



PN{x,y) + <^N{x,y)= V —D'^R{x)y^+ V D^R{x-^)y'-, 

^ — ' Q:! — ^ ot\ 

|Q;|<Ar-l \(x\=N 

where ^ is a point on the line segment from to y. Two estimates that will 
be used in the proof are 

^This is because it is differentiable away from the singularity, and its gradient is 
bounded. 
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Dl^Rix) < CII^Eir'^-l^l (5) 

for X ^ and 

\R{x-y)-PN{x,y)\<C\\x\\-<^''+'''>\\y\\^ (6) 

for \\y\\ < \\x\\/2. The estimate ([5]) can be verified by differentiating the 
Riesz kernel R{x) = 2;j/||a;||'^'*'^, and ([6]) follows by applying (j5]) to the 
remainder in Taylor's theorem. 

Theorem 3.2. If f ^ Lj^fCRf^), then the component transforms of f 
TZifix) = Cd lim / f{y)^N{x - y)dy 

^ J \\x-y\\>e' 

satisfy the bound 

|7^i/(a;)| < C7(l + 
for some < 5 < 1 — e. 

Proof. The upper bound was proved in Proposition 13.11 so it suffices to 
verify that TZif{x) has sufficient decay for large. In order to prove the 
bound, we divide the domain of the integral defining TZif into the regions 



D^ = {\\y\\ < \\x\\/2}, 

D2 = {\\y-x\\ > ||a;||/2,||y|l > ||a;|l/2}, 

Ds = {\\y-x\\ < \\x\\/2}. 

The regions Di and are disjoint open balls centered at the origin and 
X respectively, and D2 is the complement of their union. Since we assume 
||a;|| > 0, these sets form a partition of R*^. First, we consider the region Di. 
Using equation ([6]) we get 

/ \f{y)<^N{x,y)\dy<C [ \fiy)\\\yt\\x\\-''-''dy, 

JDi JDi 

and since 2 < ||a;||/||y||, the inequality 



\f{y)<^N{x,y)\dy < C\\x 



\-d-N+e+S 



\N-e-S 



/(•) 



8 



holds for any positive 6. Furthermore, due to the decay of /, the norm in 
the right hand side of the above equation is finite. 

Next, consider the region D2, where we may apply equation ([5]) to obtain 



y||''||a;|| dy 



/ \f{y)^N{x,y)\dy<C I (^||x|r^ + V 

and note that for k < N — 1, we have N — {e-\-6) — k > 0, since < e + 5 < 1. 
Additionally, 2||y|| > ||a;|| in this region, so we multiply each term of the 
sum by (2||y ||/||a;||)^^('^+^)^'^ to obtain the required bound: 



\f{y)'^N{x,y)\dy<C\\x\\-'' / \f{y)\ 
D2 Jd2 ^ 11^ 



y\ 



N-e-5 



dy 



< C\\x 



\-d-N+<E+S 



/(•) 



Finally, the region D3 contains the singularity of the Riesz kernel, and 
we further divide the integral: 



lim / /(y)$Ar(a;,y)dy = lim / f{y)R{x - y)dy 

f{y)PNix,y)dy, 

D3 

where B{x,e') denotes the ball of radius e' centered at x. The second part 
is bounded using techniques similar to those used in the region D2: 



N-l 

/ \f{y)PN{x,y)\dy<C \f{y)\Y^\\yf 



\x\\~'^~''dy 



<C [ \f{y)\\\x\ 
< C\\x\\-^-'^+'+^ 



d \\y\ 



N-l 



X 

N-e-S 



dy 



/(• 
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The first part is bounded using the mean value theorem along with the 
estimate for the decay of the derivatives of /: 



/ fiy)Rix-y)dy = [ ^-^^^^^^P^dy 

lD3\B(x,e') JDz\B(x,e') 11^^ " ^11 \\x - yT 

<C sup ||V/(t)|| / Wx-yW-'^+^dy 

teB{x,\\x\\/2) J D-i 

< C{l + \\x\\)-'^-^+\ 

As this bound is independent of e', the result follows. □ 

As previously stated, the Riesz transform preserves both decay and van- 
ishing moments, and both are vital for deriving decay estimates for the 
higher-order Riesz transforms. The former property was the subject of the 
last theorem, and we now address the latter; i.e., we will show that for any 
/ G Liv(R'^), ^if has vanishing moments up to degree A^ — 1. Basically, the 
idea is to consider both the Fourier and spatial domain formulations concur- 
rently. To wit, vanishing moments can be interpreted both as the order of 
the zero of a function at the origin in the Fourier domain and as an integral. 

Proposition 3.3. For f G Ln(R'^) and |/3| < A^, we have 

=o(||cj||^-i-|^l) as u^O. 

Proof. The decay estimate on / implies that {■)^f is in L^CHV^) for < N, 
so / has continuous derivatives up to order — 1 and 

D^f = C[i-ff]\ 

cf. [20^ Theorem 1.1.7]. The right-hand side can be expressed as the integral 
Df^f{u) = C j x^f{x)e^''-'^dx, 

and the vanishing moment condition on / implies that D^f{0) = for 
|/3| < A'^. Therefore, Taylor's theorem implies that 



L — ^ re 



a! 

\ol\=N 



for some ^ on the line segment from to cj. The result now follows from 
the fact that ^ as | ^ 0. □ 
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Theorem 3.4. Given any f € LnCR!^), TZif has vanishing moments up to 
degree N — 1. 

Proof. The proof rehes on considering two representations oiTZif. First we 
have 

TZTfiu) = C j nif{x)e^''-'^dx, 

so 

Dl^nTfiuj) = C j xl^^Z^f{x)e^''■'^dx. 



We also know that 



and hence 



\\UJ\\ 



D^^nTfio) = 

for |/3| < by Proposition 13.31 Putting the two together, we have 



= Df^nTfio) = c j x^n^f{x)dx 



□ 



Let us conclude this section by considering some examples. To see how 
our results can be applied, we will first examine the Hermitian wavelets in 
one dimension, which we recall are derivatives of the Gaussian, 

where C„ is a normalization constant. Clearly, these wavelets are infinitely 
smooth and decay faster than any polynomial. Therefore, we must look 
at their vanishing moments to determine the spaces L^r in which they lie. 
Being derivatives of the Gaussian, integration by parts shows that V'n has 
n vanishing moments, so ^„ £ Ln{R). Applying the results of this sec- 
tion, we conclude that the Riesz transform of ipn, denoted by T^Vn, satisfies: 
\TZipn{x)\ < C(l + Ixl)"" and T^V'n has n vanishing moments. Note that 
since we are only interested in the essential decay of the Riesz transform, 
we have bounded the term jxl^^*^ by \x\ in our estimate. However, the more 
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precise decay rate stated in Theorem 13. 21 will be important for showing how 
decay is preserved for higher-order Riesz transforms. 

As a second example, consider the two dimensional Poisson kernel 

p{x;s) 



27r(s2 + ||a;P)3/2' 



which was studied by Felsberg and Larkin [14| [T^. The components of the 
Riesz transform of this function are computed to be 

TZipix; s) = C, ^' 



||a;||2)3/2^ 



which have decay comparable to ||ic||^2 direction Xi. Here, we can 

again see that while the Poisson kernel is smooth, its Riesz transform has 
slower decay due to the fact that p{x; s) has no vanishing moments. 



4 Higher-order Riesz transforms 

In the previous section, we showed that a component transform of a func- 
tion in Ljv(R'^) has a rate of decay similar to the original function and a 
comparable number of vanishing moments. With this in mind, we would 
like to extend these results to higher-order Riesz transforms by applying 
the results for the single order transform in an iterated fashion. In order 
to accomplish this, we will need to impose conditions which are stronger 
than those required to be in L]\f{'R.'^). To make this precise, we recall the 
definition of the Sobolev spaces on R"^. 

Definition 4.1. Let F*(R'^) denote the Sobolev space of order s > on 
R"^, i.e. 

H'(R^) = {/ G L\li^) : f{-){l + 11-11^)^/2 e L2(R'^)} . 

The additional properties that we required of the original wavelet are 
given in terms its derivatives. In particular, the derivatives of the wavelet 
should also have decay and vanishing moments. 

Definition 4.2. For TV G N and M G N we define the space Ln,m{'R'^) to 
be the collection of / : R"^ R satisfying: 

1. / G H'(R'^) n C(R'^) for some s> M + d/2 

2. There are constants C > and < e < 1 such that for any < |q;| < 
M we have 
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(a) \D'^f{x)\ < C{1 + for < [q[ < M 

(b) / x'^D°'f{x)dx = for < |q;| < M - 1 and \P\ < N + \cx\ 

Notice that the first property in the definition implies -L7v,m(R'^) is con- 
tained in C^(R'^) by the Sobolev embedding theorem. Our definition also 
guarantees that Lj\[^]\j(R'^) is a subset of LAr(R'^) for M > 1, and we will 
show that the component transforms map Liv,M(R''^) to Ln^m-i(R'^)- Hence 
the iterated transforms of a function in L7v,Af(R-'^) will have a rate of decay 
comparable to that of the the single transform of a function in Ljv(R'^). 
The proof relies on the iterated application of Theorem 13.21 and Theorem 
13. 4t however, we must first establish that the transforms commute with 
differentiation. 

Lemma 4.3. If f e LnmC^"^) for some M > 1, then 7^^L'"/ = D^TZi f 
for any < \a\ < M — I. 

Proof. First, since / € for some s > M + d/2, the Sobolev embedding 
theorem implies / G C^^(R'^). Additionally, 

SO TZif is also in H"^ and consequently in C^^(R'^). Therefore for any fixed 
< \a.\ < M — 1, TZiW^f and D^TZif are continuous functions. Further- 
more, they are equivalent as distributions since they have the same Fourier 
transform uji/\\uj\\uj"f. Therefore they are equivalent as functions, [11', The- 
orem, p. 27]. □ 

Combining this lemma with the results of the previous section, we can 
see how the defining properties of L^^m ensure decay of the higher-order 
Riesz transforms, as we explicitly show in the next few results. 

Lemma 4.4. If f e Ln,m(R'^) for some M > I, then TZif & Ln,m-i(R''^)- 

Proof. By Lemma 14.31 it suffices to verify the second property in the defini- 
tion of LjVjM-ilR*^)- To begin, fix a multi-index a with < \a\ < M — 1. 
Since / G Ljv^a/(R'^), we know 

1. |L>"/(a;)| < C(l + 

2. |L>"+T/(a;)| < C{1 + for |7| = 1 

3. D°^f has vanishing moments up to degree — 1 + |a|. 
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Theorem 13.21 then imphes that 



D'^lZifix) = n,D'^f{x) < C{1 + 



■d-N-\a\+e+5 



for some < S < 1 — e. Additionally, Theorem 13.41 implies that D"-TZif{x) 



Theorem 4.5. If f € Lj^^m(^'^), then for some < e < 1 and any 1 < 

\a\ < M we have TZ°'f{x) < C{1 + ||a;||)~'^~'^"'"'^ and Tl°^f{x) has vanishing 
moments up to degree N — 1. 

To illustrate this result, we consider a two-dimensional analog of the 
Hermitian wavelets: 



where /3i = /32 = 2n for some positive integer n. As before, the number of 
vanishing moments is the limiting factor when classifying ipn in terms of the 
spaces Lj\f^M- Considering their Fourier transforms, we can see that ipn is 
in Ln,n- Therefore, Theorem 14.51 implies that the nth-order Riesz transform 
of Tpn, denoted by TV^ipn, decays at least as fast as (1 -|- |[a;||)~". Moreover, 
notice that if had infinitely many vanishing moments, this theorem would 
guarantee a rate of decay faster than any polynomial. In fact, the Riesz 
transform preserves the collection of Schwartz class functions with infinitely 
many vanishing moments, as is stated formally in the following corollary. 

Corollary 4.6. If ip is a Schwartz class function and all of its moments 
vanish, then TZiijj satisfies the same properties. 

5 Application to Simoncelli wavelets 

As we have seen, spatial decay of the Riesz transform is guaranteed by 
assuming smoothness, vanishing moments, and decay of the original func- 
tion. There are several examples of wavelets that are bandlimited and whose 
Fourier transform is zero in a neighborhood of the origin (cf. [22] )• The one 
property that many wavelet constructions lack is fast spatial decay. This is 
a result of constructing wavelets in the Fourier domain using a non-smooth 
cut-off function. One way to see this is to consider a Fourier integral repre- 
sentation of a wavelet ip: 



has vanishing moments up to degree A'' — 1 -|- |q:|. 



□ 



= CnDf^e 



2/2 
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If ip has n continuous derivatives and compact support, then we can integrate 
by parts to get the decay estimate 



|^(x)| ^ 



(27r)i/2 |xp 

However, a discontinuity of a lower order derivative of ip can produce terms 
with slower decay. This estimate also shows how the size of the derivatives 
of V' play a role in the decay of V'- To give a specific example and show how 
decay can be improved, we consider the Simoncelli wavelets [IS], defined 
radially in the Fourier domain by 

1^ 0, otherwise 

Note that this wavelet has vanishing moments of all degrees and finite 
Sobolev norm of all orders. However, the lack of smoothness in the Fourier 
domain causes ip to have slow decay. With Theorem 14.51 as a guideline, we 
propose a modification to improve the decay of its higher Riesz transforms. 

If we write V' as cos(^-)) o log2(^-), where o denotes function 

composition, it is clear that we can create a smooth analog if we replace the 
indicator function x by a smooth Meyer window. In particular, let be a 
function satisfying: 

1. is identically equal to 1 on the interval [—1 + e, 1 — e] 

2. 6^ is equal to outside [—1 — e, 1 + e] 

3. has at least r continuous derivatives for some r G N 

4. e,{xf + e,{x - 2)2 = 1 for < X < 2. 

Now, if we merely replace the indicator function with 0^, we will have in- 
creased the support of of our wavelet in the Fourier domain and the basis 
it generates might no longer be a tight frame. To compensate, we rescale 
so that support is maintained within [— vr, vr]'^. Let us denote the modified 
wavelet by ip[, which in the Fourier domain is defined radially by 



(0,cosg.))olog2(^-) 



The decay of this new function and its transforms can be determined by the 
results of the previous sections. 
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Table 1: Examples of finite smootiiness Meyer windows: Replacing 
G by Gn in equation ^ and using the resulting Hn,e in equation ([8]), we 
obtain a 0„ ^ in C"(R). 



n 



Gn{x) 



3 

4 
5 



357L(_1^7 I 3^5 
64 V 7-^ "T 5-^ 



3157r /I 9 
512 Vg-^ 



4 7 , 6 5 



35' 

gX^ + X + 3y|) 



69371/ ]_ 11 I 5^9 

1024 V 11-^ "T 9-^ 



5^-3 -L T- 256 ^ 
3 693' 



Proposition 5.1. The value of r defining the smoothness of Of, will deter- 
mine the decay of ip',,, and the higher Riesz transforms of tp[ will maintain a 
similar rate of decay. 

Here, we should point out that one of the nice properties of the Simoncelli 
wavelets is that they generate a tight frame, and importantly, our modifica- 
tion does not nullify this property. 

Proposition 5.2. The wavelet basis constructed from the modified Simon- 
celli wavelet ip',, is a tight frame. 

Proof. See Appendix [Al □ 

We conclude now with a discussion about the construction of the Meyer 
windows 0^. The starting point is a smooth, non-negative bump function g 
that is supported in [—1,1] and whose definite integral is 7r/2. From this 
function, we define 



G{u) := / g{t)dt 

J —oo 

and 

^,.,:.G(^)-f.c(^) (7) 

The Meyer window can then be constructed as 

e,{co) = cos{H,{lo)). (8) 
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Indeed, one can verify that for < a; < 2 we have 

e,{ujf + e,{uj-2f = cos^ {H,{uj)) + cos^{H,{uj - 2)) 

= cos2(F,(w)) + cos'^ iH,{uj) - tt/2) 
= cos'^ {H,{uj)) + sin2(i/,(w)) 
= 1 

In the case of finite smoothness, examples can be found by setting gn{x) = 
C(l — x^)" on [—1,1] and otherwise. For some particular examples, see 
Table [H In fact, the functions in the table are also derived in [lOj; however, 
the authors' construction is different from ours. Notice that each such Gn 
has n continuous derivatives, so the corresponding On^e and e ^i^^ ^^'^ 
have n continuous derivatives, which translates into decay of the wavelet in 
the spatial domain. 

In order to construct an infinitely smooth 0, we take another approach, 
as explained in jlOJ. Instead of starting with a bump function, we directly 
define the function G used in ([7]). Specifically, a function G with the required 
properties is given by 



2 A(a; + 1) + A(l - cj) 
where 



A(^) 



0, otherwise, 



and the result of this construction is a Meyer window with infinite smooth- 
ness. 



6 Conclusion 

In this paper, we have derived decay estimates for higher-order Riesz trans- 
forms of wavelets, where the conditions are given in terms of smoothness, 
decay and vanishing moments. In the case of a single Riesz transform, it was 
shown that the decay rate of the original function is essentially maintained 
by the component transforms, and this result is likely to be optimal. These 
results instill confidence in recently proposed steer able wavelet constructions, 
as the Riesz transform plays a large role in the design. Additionally, they 
point on the importance of specifying primal wavelets with a large num- 
ber of vanishing moments, a fast rate of decay, and a sufficient amount of 
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smoothness. In particular, our results indicate that having many vanishing 
moments is not sufficient to guarantee decay; a matched degree of smooth- 
ness in both the space and frequency domains is also required. Optimal 
conditions for the decay estimates of the higher-order transforms remains 
an open question. By using the single-order results in an iterated fashion, we 
may be imposing conditions that are stronger than necessary, and taking a 
more direct approach may provide improved results. Given the importance 
of the Riesz transform in steerable wavelet constructions, we view this as a 
topic that warrants further investigation. 



A Proof of Proposition 5.2 

Due to the support of it suffices to verify 

fcez 

for a; > 0. Using the definition, we have 

2 



CO 



and setting z = log2(2^+^a;/7r) gives 



i:|(.,c„.(|.))o,og,(^2 

feez ^ 



U! 



fcez 

This sum is evaluated by considering two cases that depend on z. First, for 



-e < 2; < e, we have 



E|V^:(2Mf = E \eeiz + k) cos [^{z + k)) 

fcez k=-i 



= cos' ^-z) + cos' [-{z - 1) j 9,{z - ly 

+ C0S2 {^{z + 1)) e,{z + lf 

+ {e,{z - If + 9,{z + \f) cos' - 1)) 



= sm' 



1. 
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Next, for e < 2; < 1 — e, we have 
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